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Localization length of nearly periodic layered metamaterials 
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We have analyzed numerically the localization length of light £ for nearly periodic arrangements of homo¬ 
geneous stacks (formed exclusively by right-handed materials) and mixed stacks (with alternating right and 
left-handed metamaterials). Layers with index of refraction m and thickness Li alternate with layers of index 
of refraction ri2 and thickness L2. Positional disorder has been considered by shifting randomly the positions 
of the layer boundaries with respect to periodic values. For homogeneous stacks, we have shown that the lo¬ 
calization length is modulated by the corresponding bands and that £ is enhanced at the center of each allowed 
band. In the limit of long-wavelengths A, the parabolic behavior previously found in purely disordered systems 
is recovered, whereas for A -C Li + L2 a saturation is reached. In the case of nearly periodic mixed stacks 
with the condition \niL\\ = |n2-I/2|. instead of bands there is a periodic arrangement of Lorenztian resonances, 
which again reflects itself in the behavior of the localization length. For wavelengths of several orders of magni¬ 
tude greater than L\ + L2, the localization length £ depends linearly on A with a slope inversely proportional to 
the modulus of the reflection amplitude between alternating layers. When the condition |niLi| = |n2-L2| is no 
longer satisfied, the transmission spectrum is very irregular and this considerably affects the localization length. 


I. INTRODUCTION 

During the last decades, a new type of artificial materials, 
the so-called left-handed metamaterials (LH), have attracted a 
great deal of attention. They present negative indices of re¬ 
fraction for some wavelengths i] , with considerable applica¬ 
tions in modern optics and microelectronics M. Metama¬ 
terials can resolve images beyond the diffraction limit [6, 7], 
act as an electromagnetic cloak iIb UToIi . enhance the quantum 
interference 111] or yield to slow light propagation iflZl . 

Regarding the localization length in disordered systems, the 
presence of negative refraction in one-dimensional (ID) dis¬ 
ordered metamaterials strongly suppresses Anderson localiza¬ 
tion 0. As a consequence, an unusual behavior of the lo¬ 
calization length £ at long-wavelengths A has been observed. 
Asatryan el al. reported a sixth power dependence of £ with 
A under refractive-index disorder 00 instead of the well- 
known quadratic asymptotic behavior £ ~ A 2 fl6l - fl9ll . Re¬ 
cently, Mogilevtsev et al. 120 11 have also found a suppres¬ 
sion of Anderson localization of light in ID disordered meta¬ 
materials combining oblique incidence and dispersion while 
Torres-Herrera et al. [21] have developed a fourth order per¬ 
turbation theory to resolve the problem of non-conventional 
Anderson localization in bilayered periodic-on-average struc¬ 
tures. The effects of polarization and oblique incidence on 
light propagation in disordered metamaterials were also stud¬ 
ied in Ref. B22ll . 

In this article, we calculate numerically the localization 
length of light £ for a one-dimensional arrangement of layers 
with index of refraction n\ and thickness L \ alternating with 
layers of index of refraction 77.2 and thickness L 2 . In order to 
introduce disorder in our system, we change the position of 
the layer boundaries with respect to the periodic values main¬ 
taining the same values of the refraction indices n\ and ri 2 . 
This is the case of positional disorder, in contrast to the com¬ 
positional disorder where there exist fluctuations of the index 
of refraction ll23ll . 

Two structures will be analyzed in detail: homogeneous 
stacks (H), composed entirely by the traditional right-handed 
materials (RH) with positive indices of refraction, and mixed 


stacks (M) with alternating layers of left- and right- handed 
materials. For the sake of simplicity, the optical path in both 
layers will be the same, that is, the condition |mLi| = \n 2 L 2 \ 
is satisfied in most of the work. These periodic-on-average 
bilayered photonic systems have already been studied analyt¬ 
ically by Izrailev et al. lf24l [25^ . These authors have devel¬ 
oped a perturbative theory up to second order in the disorder to 
derive an analytical expression for the localization length for 
both H and M stacks. In our case, we have obtained two equa¬ 
tions for the localization length £ as a function of the wave¬ 
length A from our numerical results. For H stacks, a quadratic 
dependence of £ for long-wavelengths is found, as previously 
reported in the literature. On the other hand, the localization 
length saturates for lower values of A. An exhaustive study 
of £ in the allowed and forbidden bands (gaps) of weakly dis¬ 
ordered systems will be carried out. We will show that the 
localization length is modulated by the corresponding bands 
and this modulation decreases as the disorder increases. For 
low-disordered M stacks and wavelengths of several orders 
of magnitude greater than the grating period A = L\ + L 2 , 
the localization length £ depends linearly on A with a slope 
inversely proportional to the modulus of the reflection ampli¬ 
tude between alternating layers. 

The plan of the work is as follows. In Sec. II we carry out an 
exhaustive description of our one-dimensional disordered sys¬ 
tem and the numerical method used in our localization length 
calculations. A detailed analysis of £ in the allowed bands and 
gaps of homogeneous stacks is performed in Sec. Ill where a 
practical expression for the localization length as a function 
of A and the disorder is derived. In Sec. IV we calculate £ 
for mixed stacks of alternating LH and RH layers. A linear 
dependence of the localization length at long-wavelengths is 
found for low-disordered M stacks. Finally, we summarize 
our results in Sec. V. 


II. SYSTEM DESCRIPTION AND NUMERICAL MODEL 

Let us consider a one-dimensional arrangement of layers 
with index of refraction ni alternating with layers of index of 
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N- 1 layers 



where the index j runs from 1 to N and the coefficients A } 
and Bj can be written as 
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Figure 1. A periodic arrangement of layers with index of refraction 
m and thickness Li alternating with layers of index of refraction ri 2 
and thickness L 2 . The grating period is A = Li + Z/ 2 - 


refraction ri 2 . The width of each one is the sum of a fixed 
length Li for i = 1,2 and a random contribution of zero mean 
and a given amplitude. The wave-numbers in layers of both 
types are ki = um,i/c, where lo is the frequency and c the 
vacuum speed of light. As previously mentioned, the grat¬ 
ing period of our system A is defined as the sum of the aver¬ 
age thicknesses L\ and L 2 of the two types of layers, that is, 
A = Li + L' 2 - We have introduced the optical path condition 
|niLi| = I 77 . 2 A 2 I for simplicity (in the case of left-handed lay¬ 
ers rii < 0 , so the absolute value has been written to consider 
these type of materials). Without disorder, each layer would 
be limited by two boundaries x^ and where N is the 
total number of boundaries. The periodic part of the system 
considered is schematically represented in Fig. [I] 

In the presence of disorder, the position of the correspond¬ 
ing boundaries are 


The parameters rj-ij, which are the reflection amplitudes be¬ 
tween media j — 1 and j, are given by 


r 7 -l ,.7 — r i,. 7 -l — 


Zj -1 — Zj 

Zj — 1 + Zj ’ 


( 6 ) 


where Zj corresponds to the impedance of layer j and can 
be be expressed for normal incidence in terms of its dielectric 
permittivity Ej and magnetic permeability fij 


Zj = 



( 7 ) 


The quantity \j-i,j entering Eqs. ((4} and ([5} is a phase 
term [261 


Aj-ij = Ayj-i = exp [ 2 ikj-i\xj — £7+11] • (8) 

Here k 3 _ 1 is the wave-number in a layer with boundaries x :j 
and Xj + i. This recurrence relation facilitates the numerical 
computation of the determinant. The initial conditions are the 
following 


x j =xf' ) +£ j 5, j = 2, ...,7V- 1, (1) 


M = 1; D 0 = 1; D_i = 0. (9) 


except for the first and the last boundary, so as to maintain 
the same total length L. The parameters £7 are zero-mean 
independent random numbers within the interval [—0.5, 0.5]. 
Throughout all our calculations, we have chosen values of the 
disorder parameter <5 less than L \ and L 2 . 

For each L, we calculate the transmission coefficient of our 
structure T and average its logarithm. In T, over 800 disorder 
configurations. Then, we obtain numerically the localization 
length £ via a linear regression of In T (23] 


lim 

L —>00 


(InT) 
2 L 


1 

l 


( 2 ) 


Here, the angular brackets (...) stand for averaging over the 
disorder. We choose 6 values of the total length L to perform 
the linear regression of Eq. Q. The localization length £ is 
evaluated as a function of the disorder parameter 6 and the 
frequency of the incident photon oj. 

We calculate the transmission coefficient of our system via 
the characteristic determinant method, firstly introduced by 
Aronov et al. [26]. This is an exact and non perturbative 
method that provides the information contained in the Green 
function of the whole system. In our case, the characteristic 
determinant Dj can be written as |26] 


77.7 — AjDj -1 BjDj-2, ( 3 ) 


The transmission coefficient of our structure T is given in 
terms of the determinant Dm by 

T=\D n \~ 2 . (10) 

III. LOCALIZATION LENGTH FOR HOMOGENEOUS 
STACKS 

Before dealing with mixed stacks, we present results for 
low-disordered homogeneous systems with underlying peri¬ 
odicity, which has not been previously studied. In this sec¬ 
tion we perform a detailed analysis of the localization length 
£ in the allowed bands and in the forbidden gaps of disor¬ 
dered H stacks as a function of the disorder <5, the incident 
wavelength A and the reflection coefficient between alternat¬ 
ing layers |rj_ij| 2 . 

As it is well known, in the absence of disorder the trans¬ 
mission spectrum of right-handed systems presents allowed 
and forbidden bands whose position can be easily determined 
via the following dispersion relation obtained from the Bloch- 
Floquet theorem |27[| 

cos(/3A) = cos(/tiLi) cos(fc 2 A 2 ) 

- | sin(fciLi)sm(fc 2 L2), ( !1 ) 
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Figure 2. (a) The transmission coefficient T and (b) the parameter 
cos(/3A) versus the frequency cj for the homogeneous periodic sys¬ 
tem described in the text (99 layers). 


where f3 is the Block wave-vector. When the modulus of the 
right-hand side of Eq. O is greater than 1, /3 has to be taken 
as imaginary. This situation corresponds to a forbidden band. 
Taking into account the condition |niLi| = \n 2 L 2 \, Eq. (ITO 
reduces to 

cos(^A) = cos 2 (fciLi)-i sin 2 (fciTi). (12) 


is, an arrangement of layers with random boundaries and al¬ 
ternating indices of refraction n\ and n 2 . Several features are 
evident in the figure. For long-wavelengths, one observes a 
quadratic asymptotic behavior, as can be compared with the 
dotted line j 1 6h 1 911 . An in-deep numerical analysis of the co¬ 
efficient characterizing this dependence has been performed. 
To this aim, 20 different H stacks were considered and the 
following expression for the localization length was found 

A 2 

£ ~ 0.063———, for A ->• oo (14) 

Ao P r-o- 

where A op = niLi+n 2 L 2 is the optical path across one grat¬ 
ing period A. All the lengths in Eq. llT4l > are expressed in units 
of A. In the opposite limit of short A, the localization length 
£ saturates to a constant value iTl . [28h . Our numerical results 
have shown that this constant is proportional to the inverse of 
the reflection coefficient between alternating layers |r| 2 , that 
is, 

£^i for A —»0. (15) 

Izrailev et al. | |24l[25|| have developed a perturbative theory 
up to second order in the disorder to calculate analytically the 
localization length in both homogeneous and mixed stacks. 
This model is quite general and is valid for both quarter stack 
medium (mainly considered in our work) and systems with 
different optical widths. Assuming uncorrelated disorder and 
random perturbations with the same amplitude in both layers 
(the main considerations in our numerical calculations) one 
can easily derive the following analytical expression for £ at 
long-wavelengths from Izrailev’s formulation 


On the other hand, when cos(/3A) is equal to unity, the inci¬ 
dent frequency co is located at the center of the m-th allowed 
band, ! . After some algebra, we obtain from Eq. (1 1 21 


co 


( m ) — 
C 


m7T 




(13) 


, _ ^2 2A 2 

{Z\ - Z 2 ) 2 7r 2 (nf + rc 2 )<5 2 

For similar values of the layer impedances Z\ ~ Z 2 , the first 
term in Eq. ( IT6l > can be approximated by 1 /4r 2 and n 2 + ~ 

2A 2 p , so Eq. (IT6l > reduces to 


Let us first consider a periodic H stack formed by 50 lay¬ 
ers of length Li = 52.92 nm and index of refraction n\ = 
1.58 alternating with 49 layers of length L 2 = 39.38 nm and 
n 2 = 2.12. The total size of our structure is 4.57 /rm and the 
reflection coefficient between alternating layers 0.05259. Fig. 
[2]' a) represents the transmission coefficient T as a function of 
the frequency oj to illustrate its behavior. Also shown are the 
center of each allowed band calculated via Eq. (Il3l) . There are 
99 peaks in each band so they can hardly been resolved on the 
scale used. Moreover, in Fig. [2jb) the parameter cos(/3A) is 
plotted versus the frequency co for this periodic system. The 
first gap and the first allowed band have been shown for a bet¬ 
ter comprehension. 

A systematic numerical simulation of a realistic system 
with 50000 layers has been carried out. The parameters are 
the same as in the previous example. In Fig. [3] we represent 
the localization length £ versus the wavelength A for different 
values of the disorder parameter S (shown in the legend of the 
figure). The dashed line corresponds to ’’total disorder”, that 


£ 


(±-)* 

\47 T 2 J A 2 p 7’ 2 <5 2 


0.025 


A 2 

A 2 p r 2 <5 2 ’ 


for A —> oo 


(17) 


which is similar to our numerical expression Eq. d. 

The randomness only affects partially the periodicity of the 
system, which manifests in the existence of bands and gaps. 
The localization length depends on the position in the band 
and on the disorder. The modulation of £ by the bands can 
be clearly appreciated in Fig. [3] These results are consistent 
with other published works on this topic j29ll3^1 . Recently, 
Mogilevtsev et al. [29] have reported that the photonic gaps 
of the corresponding periodic structure are not completely de¬ 
stroyed by the presence of disorder while Luna-Acosta et al. 
l30ll have shown that the resonance bands survive even for rel¬ 
atively strong disorder and large number of cells. 

Having a close look into the first gap in Fig.[3j °ne observes 
that the localization length is practically independent of the 
disorder S. In order to visualize this effect. Fig. 0 represents 
(a) the first and (b) the second gaps depicted in Fig. [3] As 
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Figure 3. Localization length £ versus the wavelength A for different 
values of the disorder parameter S. The H stack corresponds to the 
arrangement represented in Fig. [2] but now 50000 layers have been 
considered. The dashed line stands for the "’total disorder” case. All 
lengths are expressed in units of the grating period A. 


Figure 5. Three-dimensional graphs of the localization length £ ver¬ 
sus the wavelength A and the disorder S for (a) the first and (b) the 
third allowed bands. The FI stack is the same as in Fig. [2] All lengths 
are expressed in units of the grating period A. 



Figure 4. Localization length £ versus the wavelength A for (a) the 
first and (b) the second gaps depicted in Fig. [3] 


mentioned, the dependence of £ with the disorder is almost 
negligible in the first gap. When the wavelength is similar to 
the grating period A, the influence of the disorder is greater, 
as can be easily deduced from simple inspection of Fig. 0b). 

Let us now focus on the allowed bands and study in detail 
the behavior of the localization length in these regions. To 
this aim, three-dimensional (3D) graphs of £ versus the wave¬ 
length A and the disorder 5 have been plotted in Fig.0for (a) 
the first and (b) the third allowed bands (see again Fig. 0. All 
this magnitudes have been normalized to the grating period A. 
The localization length £ is enhanced in a small region around 
the center of each allowed band. A similar result was found by 
Hernandez-Herrejon et al. m who obtained a resonant effect 
of £ close to the band center in the Rronig-Penney model with 
weak compositional and positional disorder. This increase in 
the localization length is due to emergence of the Fabry-Perot 
resonances associated with multiple reflections inside the lay¬ 
ers from the interfaces 112411251 1321. In particular, for homoge¬ 


neous quarter stack systems, the Fabry-Perot resonances arise 
exactly in the middle of each allowed band where /3 vanishes 
1124112511 . The saturation of £ for short-wavelengths is also ap¬ 
preciated in these 3D images. 

Up to now, H stacks with the same optical path in layers of 
both types have been considered, that is, arrangements verify¬ 
ing the condition \n\Li\ = |rz 2 -^ 2 1 in the absence of disorder. 
As a consequence, the transmission spectrum T of the corre¬ 
sponding periodic system presented a symmetric distribution 
of allowed bands and gaps (as previously shown in Fig. 0. 
What happens in the case of a non-symmetric band distribu¬ 
tion, that is, when the condition \n\Li\ = |tt. 2 -^ 2 | is n °t satis¬ 
fied? To answer this question, we have plotted the transmis¬ 
sion coefficient T (Fig. 0a)) and the parameter cos(/3A) (Fig. 
0b)) versus the frequency oj for a periodic H stack formed by 
50 layers of length L\ = 52.92 nm and index of refraction 
ni = 1.58 alternating with 49 layers of length L 2 = 28.80 nm 
and ri 2 = 2.12. Note that the condition \n\Li\ = |ri 2 -^ 2 1 is 
no longer held, so the band structure is asymmetric. Accord¬ 
ingly, the localization length £ shown in Fig. 0c) presents an 
irregular form in the allowed and forbidden bands. As in the 
symmetric case, no band modulation exists for high disorders 
and the quadratic asymptotic behavior for long-wavelengths 
is also verified. Moreover, the peaks in the localization length 
due to Fabry-Perot resonances still can be appreciated, al¬ 
though they are no longer in the center of the bands ll24l [25tl . 
A total number of 50000 layers was considered in our local¬ 
ization length calculations. 


IV. LOCALIZATION LENGTH FOR MIXED STACKS 

Once analyzed in detail the behavior of the localization 
length £ for homogeneous systems, let us now deal with M 
stacks composed of alternating LH and RH layers. 

In our numerical calculations we have considered a peri¬ 
odic M stack formed by 50 layers of length L\ = 52.92 nm 
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Figure 6. (a) The transmission coefficient T and (b) the parameter 
cos(/3A) versus the frequency oj for the asymmetric periodic H stack 
described in the main text (99 layers) and (c) the corresponding lo¬ 
calization length £ versus the wavelength A for different disorder pa¬ 
rameters 5 (50000 layers). 
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Figure 7. (a) The transmission coefficient T and (b) the parameter 
cos(/3A) versus the frequency c j for the mixed periodic system de¬ 
scribed in the text (99 layers). 


and index of refraction n\ = -1.58 alternating with 49 layers 
of length L 2 = 39.38 nm and n 2 = 2.12. Again, the condition 
\niLi\ = \n 2 L 2 1 has been imposed. Note that this arrange¬ 
ment has similar parameters than the one depicted in Sec. Ill, 
but now n\ is negative. This change of sign results in a se¬ 
vere modification of the transmission coefficient T, as we will 
show immediately. For this the periodic system, Fig.[7]repre- 
sents (a) the transmission coefficient T and (b) the parameter 
cos(/3A) versus the frequency u of the incident light. Unlike 
the H stack case, no allowed bands exist and practically the 
entire transmission spectrum is formed by gaps. A set of pe¬ 
riodically distributed Lorentzian resonances is found instead. 
The position of the center of each resonance is given by Eq. 
o. that is, the center of the allowed bands in homogeneous 
systems. 

In respect to the localization length, positional disorder was 
introduced as explained in Sec. II. As previously considered, 
the total number of layers in our numerical calculations was 
50000 and the number of disordered configurations to average 
the logarithm of the transmission coefficient was 800. The re¬ 
sult is shown in Fig. [8] where the localization length £ is rep¬ 
resented versus the wavelength A for different values of the 
disorder parameter <5. The dashed line corresponds to the ’’to¬ 
tal disorder” case. Again, for long-wavelengths a quadratic 
asymptotic behavior of £ is found, but now a region where the 
localization length is proportional to A exists. We will turn to 
this point in the next figure to quantify the slope of this lin¬ 
ear dependence. As it is noticed, the Lorentzian resonances 
associated with multiple reflections in the layers modulate the 
shape of £ and this modulation decreases as the disorder in¬ 
creases. Moreover, the saturation of the localization length for 
low-wavelengths can also be appreciated. As in the H stack 
case, the constant where £ saturates is proportional to the in¬ 
verse of the reflection coefficient between alternating layers 



Figure 8. Localization length £ versus the wavelength A for different 
disorder parameters 8. The M stack corresponds to the one repre¬ 
sented in Fig. [7] but here 50000 layers have been considered. The 
dashed line stands for the ’’total disorder” case. 


The linear dependence of £ with the wavelength A has been 
exhaustively studied by our group to find a simple analytical 
expression for the localization length in this region. More than 
30 different M stacks have been simulated and we have arrived 
at the following empirical equation 


a 


A 

6 A op |r| 


= aX, 


( 18 ) 


where £, A and A op are expressed in units of the grating pe¬ 
riod A. In Fig. [9] our numerical calculations of the slope a 
versus |r| have been plotted for several values of A op , trian¬ 
gles (1.25), squares (3.25) and circles (7.55). The solid lines 
correspond to the results obtained via Eq. dH. One notices a 
good degree of validity for a wide range of |r| values. 

Finally, let us now consider an asymmetrical M stack where 
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Figure 9. Numerical calculations of the slope a versus |r| for several 
values of A op (expressed in units of the grating period A). The solid 
lines correspond to the results obtained via Eq. GD 



(o (THz) x (units of A) 


Figure 10. (a) The transmission coefficient T and (b) the parame¬ 
ter cos(/3A) versus the frequency uj for the asymmetric periodic M 
stack described in the main text (99 layers) and (c) the corresponding 
localization length £ versus the wavelength A for different disorder 
parameters <5 (50000 layers). 


the condition |niLi| = \n 2 L 2 \ is no longer satisfied. In Fig. 
flOlwe have represented (a) the transmission coefficient T and 
(b) the parameter cos(/3A) versus the frequency uj for a peri¬ 
odic M stack formed by 50 layers of length of length L\ = 
52.92 nm and index of refraction ni = -1.58 alternating with 
49 layers of length L 2 = 28.80 nm and n 2 = 2.12. Note the 
strong difference between this transmission spectrum and the 
symmetrical one (see Fig. [TJa)) where a set of periodically 
distributed Lorenztian resonances exists. Despite this fact, the 
localization length £ shown in lug. 11 Ol e) presents a region of 
linear dependence with the wavelength, as in the symmetric 
case. However, Eq. (fl8l) cannot be used to evaluate the local¬ 
ization length in this region. 


V. DISCUSSION AND CONCLUSIONS 


We have analyzed numerically the localization length of 
light £ for homogeneous and mixed stacks of layers with index 
of refraction ±|m| and thickness L\ alternating with layers 
of index of refraction \ri 2 \ and thickness L 2 . The positions of 
the layer boundaries have been randomly shifted with respect 
to ordered periodic values. The refraction indices n\ and JI 2 
present no disorder. 

For H stacks, the parabolic behavior of the localization 
length in the limit of long-wavelengths, previously found in 
purely disordered systems 03-B, has been recovered. On 
the other hand, the localization length £ saturates for very low 
values of A. The transmission bands modulate the localiza¬ 
tion length £ and this modulation decreases with increasing 
disorder. Moreover, the localization length is practically inde¬ 
pendent of the disorder S at the first gap, that is, it has a very 
low tendency in this region. We have also characterized £ in 
terms of the reflection coefficient of alternating layers |r| 2 and 
the optical path across one grating period A op . Eq. ( 1 1 4| i has 
been proved to be valid for a wide range of |r| 2 values, that is, 
from transparent to opaque H stacks. It has also been shown 
(see Fig. [5]) that the localization length £ is enhanced at the 
center of each allowed band. 

When left-handed metamaterials are introduced in our sys¬ 
tem, the localization length behavior presents some differ¬ 
ences with respect to the traditional stacks, formed exclusively 
by right-handed materials. For low-disordered M stacks and 
wavelengths of several orders of magnitude greater than the 
grating period A, the localization length £ depends linearly on 
A with a slope inversely proportional to the modulus of the 
reflection amplitude between alternating layers |r| (see Eq. 
(fl8l) ). As in the H case, £ saturates for low-wavelengths, being 
this saturation constant proportional to the inverse of |r| 2 . 

If we take into account losses, there is an absorption term 
whose absorption length £ a b s is El 


£abs 


A 

2na ’ 


(19) 


where a is an absorption coefficient. The inverse of the total 
decay length is the sum of the inverse of the localization length 
£ plus the inverse of the absorption length £ a b s . Note that 
£abs is proportional to A, so, for low-disordered M stacks and 
weak absorption metamaterials, the final expression for the 
localization length £ in the linear region can be written as 


£ 


A 

6A op |r| + 27rcr 


( 20 ) 


In the case of both homogeneous and mixed stacks with 
non-symmetric band distribution, that is, when the condition 
|niLi| = I 712 L 2 I is not satisfied, the localization length £ 
presents an irregular form in all the transmission spectrum. 
These changes in £ are more sensitives in mixed stacks than 
in homogeneous structures. 
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